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Abstract 

Rietveld profiles, Mossbauer spectra and x-ray absorption fine structure (XAFS) 
were analyzed through the structural phase transition at Ts — 355 K in HoBaCo407. 
Excess of the oxygen content over "O7" was avoided via annealing the samples 
in argon flow at 600 °C. Space groups (S.G.) Pbn2ic and PQ^mc were used to 
refine the structure parameters in the low- and high-temperature phases, respec- 
tively. Additionally, the Cmc2i symmetry was considered as a concurrent model 
of structure of the low-temperature phase. In the high-temperature phase, se- 
vere anisotropy of thermal motion of the major part of the oxygen atoms was 
observed. This anisotropic motion turns to be quenched as the sample is cooled 
below Tg. The variation of quadrupole splitting e near Ts is not similar to a 
steplike anomaly frequently seen at the charge-ordering transition. We observe 
instead a dip-like anomaly of the average e near Tg. Narrow distribution of 
the electric field gradient (EFG) over different cobalt sites is observed and ex- 
plained on the basis of point-charge model. XAFS spectra show no evidence of 
significant difference between YBaCo407 (T > Ts) and HoBaCo407 (T < Ts). 
The origin of the transition at Ts is ascribed to the condensation of the libra- 
tion phonon mode associated with the rigid-body rotational movements of the 
starlike tetrahedral units, the building blocks of kagome network. It is shown 
that the condensation of the libration mode is not compatible with translation 
symmetry for the hexagonal S.G., but compatible for the orthorhombic S.G. 
The orthorhombic lattice parameters and EFG components (Vxx, Vyy, Vzz) vary 
smoothly with temperature at approaching Tg and closely follow each other. 
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Introduction 

The cobaltites REBaCo407 have displayed a structural phase transition 
manifestative through stepwise changes of physical properties. Peak-like anoma- 
lies of specific heat were observed for RE=Lu, Yb, Tm, Er in the range of tran- 
sition temperatures Tg between 160 K (Lu) and 280 K (Er)[l|,[3|. Magnetization 
[l|, Seebek coefficient j4| and resistivity [1, |3| drop abruptly as temperature 
increases through the RE-specific value of Tg. Because of the mixed valence 
of Co (averaged valence " -1-2.25") there may occur different distributions of 
charges in the low- and high-temperature phases. Therefore, the transition was 
conjectured to be driven by the charge ordering in the Co subsystem|l|. 

REBaCo407 belongs to the family of swedenborgite, whose structure was 
solved by Pauling et alQ in 1935. The swedenborgite SbNaBe407 was de- 
scribed as a net built up of starlike Be40i3 clusters[f3i, and similar C04O13 
stars can be found in REBaCo407. Two nonequivalent sites of Co exist in 
high-temperature hexagonal phase 0, Q , but four Co sites in low-temperature 
orthorhombic phase [H. Clearly, some redistribution of charge should be associ- 
ated with changing the symmetry. It remains under debate, however, whether 
or not the transition itself is driven by the charge redistribution. 

Both phases consist of regular 1:1 stacking of kagome and triangular layers 
of C0O4 tetrahedra (Fig. 1). In each of the phases, there are 75% of Co sites 
in kagome layers and 25% of Co sites in triangular layers. Each C04O13 star 
thus involves three tetrahedra in K-layer and one in T-layer. We found relative 
arrangement of these stars unchanged, however their exact geometry somewhat 
changes through Tg. 

From analyzing the bond valence sums in both the low- and high-temperature 
phases Huq et alfsj derived the preference of the Co site in the triangular layers 
for the Co^+ ion at T < Tg, but, contrarily, for the Co^+ ion at T > Tg. On 
the other hand, it was argued recently Q that the Co sites in triangular layers 
are generally exhibiting shorter Co-0 distances than the sites in kagome layers, 
irrespectively of the transition at Tg; the Co^^^ ions should always prefer the site 
in triangular layers. 

In this work, we investigate the phase transition in IIoBaCo407 using the 
synchrotron diffractometry, Mossbauer spectroscopy, and XAFS. It was verified 
that the substitution of 2.5% of Co with Fe did not introduce a significant 
change of Tg ~ 355 K. We observed only very small changes in Mossbauer 
spectra through Tg and we argue that no charge ordering takes place at the 
phase transition. Different symmetry groups were tested to fit the phases below 
Tg (Cmc2i and Pbn2i) and above Tg (P63TOC and P31c). From the results 
of our structural study the origin of the phase transition is understood to be 
unrelated to charge redistribution. We consider instead the double- well potential 
models [l3|, in which the transition can be explained as the condensation of a 
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Figure 1: Perspective view of the quadruple unit cell in the low-temperature phase of 
ReBaCo407. Kagome and triangular layers of C0O4 tetrahedra are marked by "K" and 
"T". The quadruple unit cell of HoBaCo407 is shown with atomic coordinates refined in S.G. 
Pbn2i. 



libration phonon mode associated with the rigid-body motion of the tetrahedra 
in kagome layers around the center of the abovementioned C04O13 stars. 



1. Experimental 

The samples of HoBaCo407 and HoBaC03.9Feo.1O7 were prepared from 
H02O3 , BaC03, C03O4 and ^''Fe203 powders using standard ceramic synthesis 
technology at 1100 °C. One sample of YBaCo407 for the XAFS study was also 
prepared using Y2O3 in the same conditions. The sample of HoBaCo407 was 
made first in Japan and exactly the same synthesis was reproduced in Russia. 
Two samples showed the reproducible values of Tg and lattice parameters. Other 
samples, doped with ^^Fe, were prepared in Russia. Final annealing in flowing 
pure Ar gas was always done at 600°C. 

The measurements of x-ray diffraction patterns for Rietveld analysis were 
conducted at the " Anomalous Scattering" beamline of the VEPP-3 storage ring 
in Siberian Synchrotron Radiation Centre (SSRC). The beamline is equipped 
with Si (111) monochromator on the primary beam and Ge (111) crystal ana- 
lyzer on the diffracted beam. Three patterns of HoBaCo407 were measured at 
300 K for three wavelengths shown in Fig. 2 (a). Fourth pattern was measured at 
380 K for A =1.5421 A. Structure parameters were refined through the analy- 



sis of full-profile x-ray diffraction intensities using FULLPROF program|llLll2l|. 
Anomalous dispersion corrections for the atomic scattering factors of Ho,Ba and 
Co were introduced into the input files of FULLPROF using the Brenann and 
Cowan data from DispAnoV2 program. 
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Figure 2: Three wavelengths used in the diffraction experiments shown on the plot of theo- 
retical atomic scattering factor corrections // and /" for Co and Ho (a); characteristic area 
of the high-resolution synchrotron x-ray diffraction patterns (Rietveld plots) at A = 1.5421 A 
in the high-temperature (b) and low-temperature (c) phases of HoBaCo407. Structure of the 
low-temperature phase was refined at T = SOOii' using either S.G. Pbn2± or Cmc2i, one of 
which (P6n2i) is shown in (c). Theoretical positions of the permitted reflections are shown 
in (c) for both Pbn2i and Cmc2i groups. Indicated by the arrow reflection (105) is among 
very weak peaks extinguished for S.G. Cmc2i, but allowed for S.G. Pbn2\. 
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The measurements of x-ray diffraction patterns in the temperature range 
between 27°C and 113°C were also conducted using a MXP21 Mac Science 
diffractometer with the following operation conditions: 50° < 26 < 60° with the 
step size of 0.02°, Cu-K„ radiation (A= 1.5405 A and 1.5443 A), V = 45 kV 
and I = 100 mA. The lattice parameters were refined using Rietveld analysis. 

The Co K-edge XANES spectra were measured at the EXAFS station of 
SSRC. The beam was monochromatized with channel-cut Si (111) monochro- 
mator. The energy resolution was 0.8 eV. 

Mossbauer spectra were collected with the velocity-reversive spectrometer 
NZ-640 (Hungary), using the regime of constant acceleration for moving source 
of '^^Co embedded in a metal (Cr) matrix. The chemical shifts are given relative 
to Fe metal. The temperature of sample was maintained with the accuracy of 
±0.1 deg. The Mossbauer source was kept at room temperature. (24±1 
°C). In this setup, the drift of the source temperature had to be taken into 
account only in the temperature dependencies of the chemical shift. Mossbauer 
spectra were measured in the range of velocities from from -4 to 4 mm/s and 
stored into 1000 channels equidistributed in velocity. In each of 27 spectra 
obtained, the number of counts per channel was approximately 10^. 

2. Results 

Differential thermal analysis results have revealed the highest temperature of 
the structural transition Tg — 355 K for HoBaCo407 among all the studied previ- 
ously cobaltites of the RBaCo407 series (R=Lu,Yb,Tb,Er,Y) (cf. Refs. J.,^,ll). 
The DTA curve collected at the cooling/heating rate of 3 K/min showed the en- 
dothermic effect at the temperature Tg of 82°C for sample heating (Fig. 3). The 
latency interval, within which the high-temperature phase can be undercooled, 
evidences the 1-st order character of the transition in agreement with the gen- 
eral renormalization group theory of the phase transitions (see, for example, the 
references in the review [lOj). 

The spacing between heating and cooling curves is collapsed in Fig. 4, where 
the temperature dependence of the lattice parameters is shown. All the data of 
Fig. 4 were obtained using the overnight program, so that the average rate of 
the temperature "sweep" was slower by two orders of magnitude than that in 
Fig.3. 

2.1. Structure refinements 

The lattice parameters below and above Tg were refined using Rietveld anal- 
ysis. To refine the lattice parameters and atomic coordinates at T < Ts we 
used previously suggested structure models based on the S.G. Cmc2i 4|, and 
Pbn2i0. The substitution of 2.5% Co with Fe has no significant effect on 
the lattice parameters. In HoBaC03.gFeo.1O7 we obtained a = 6.3013(1) A, 
b=10.9552(2) A, c=10.2225(2) A. In the Fe-doped sample HoBaC03.9Feo.1O7 
we obtained at room temperature the orthorhombic cell with a ~ 6.3006(3) A, 
b=10.9573(5) A, c=10. 2275(5) A. 
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Figure 3: The curves of differential thermal analysis in HoBaCo407 measured upon heating 
and cooling at a rate of 3 deg./min. 



The temperature evolution of the lattice parameter a and the reduced lat- 
tice parameter 6/-\/3 shows that there exist below a wide temperature range 
(AT = 30 ~ 40 K), in which the orthorhombic distortion gradually increases at 
cooling. This behavior is contrasting to the abrupt jumping of the a param- 
eter, reported for YbBaCo407[3]. The c-parameter, on the opposite, exhibits 
smaller jump at Ts in our Fig. 4, as compared to its jump reported previously 
for YbBaCo407[4]. 



Table 1. Atomic coordinates for HoBaCo407 at 300 K. 



Pbn2i 


X 


y 


z 




Cmc2i 


X 


y 


z 




Ho 


0.0014(5) 





0.6672(2) 


0.6664(2) 


0.8717(2) 


0.8720(1) 


1.4(2) 


1.3(2) 


Ba 








2/3 


2/3 


1/2 


1/2 


1.4(2) 


1.4(2) 


Col 


-0.011(1) 





-0.0055(6) 


0.0067(4) 


0.9409(4) 


0.9388(4) 


0.7(2) 


0.9(2) 


Co21 


-0.001(1) 





0.1747(2) 


0.1770(2) 


0.6996(4) 


0.6985(4) 


0.6(4) 


0.5(2) 


Co22 


0.257(1) 




0.0879(4) 




0.1891(6) 




1.1(4) 




Co23** 


0.252(1) 


0.2525(4) 


0.9239(4) 


0.9177(2) 


0.6819(6) 


0.6871(4) 


1.0(4) 


1.1(3) 


01 


0.008(4) 





-0.003(2) 


0.0122(9) 


0.2555(7) 


0.2535(7) 


1.7(7) 


1.0(5) 


021 


0.783(2) 




0.2648(14) 




0.7793(10) 




0.9(7) 




022 


-0.003(5) 





0.4860(9) 


0.4784(7) 


0.2400(13) 


0.2238(7) 


0.6(7) 


0.4 


023** 


0.735(2) 


0.763(1) 


0.7512(14) 


0.7395(7) 


0.2273(12) 


0.2600(6) 


0.5(7) 


1.6 


031 


-0.053(2) 





0.1535(11) 


0.1439(6) 


0.5116(16) 


0.5077(9) 


0.4(6) 




032 


0.224(2) 




0.0983(12) 




0.0141(22) 




1.5(9) 




033** 


0.254(3) 


0.250(1) 


0.9355(11) 


0.9193(9) 


0.5064(22) 


0.5166(8) 


1.5(9) 


**** 



*Anisotropic thermal displacement factors (ATDF) were refined for 031 and 033 in S.G. 
Cmc2i. 
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Figure 4: Evolution of the lattice cell parameters a, b/y/3 and c with temperature near the 
structural phase transition. 
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Figure 5: The conventional (main panel) and differential (inset) Rietveld plots for the 
synchrotron x-ray diffraction data collected at T=300 K (below the phase transition) in 
HoBaCo407. Main panel: observed, calculated and difference intensities near Co K- 
edgc (A = 1.6134 A). Inset: the differences /obslCoK)— /obs(tloLni) and /theor 

(CoK)- 

Ahcor(HoLjjj), where the profile Jtheor(HoLni) was transformed to the /theor(CoK) profile 
conditions. 
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**Position multiplicity 8 in the S.G. Cmc21 . All other positions in both Pbn2i and 
Cmdli groups have the multiplicity 4. 

***Refined ATDF for 031: Uii= 0.13(2), C/22=0.016(8), C/33=0.004(9), C/23=-0.006(8) 
A^, The values of and J723 a^re fixed to by symmetry. Uy are smaller than By by the 
factor of 87r'^. 

****Refined ATDF for 033: [/ii=0.033(17), C/22=0.09(3), ?733 =0.028{17), C/i2=0.01(2), 
t7l3=-0.01(2), J723=-0.05{3) A^. 

Three patterns measured on synchrotron radiation at 300 K were fitted all 
together, varying at once the different profile parameters particular of each 
pattern and the atomic coordinates common for all patterns. Results shown 
in Table 1 were obtained with the S.G.'s Pbn2i and Cmc2i. Corresponding 
rehability factors are = 3.83, i?Ll"g' = 7-4%, i^Brafg = 7-7%, i?Llfg = 9-3% 
for P6n2i, and - 3.91, i?™ = 7.1%, R},^^^^ = 7.2%,R°,f^^^ = 9.9% for 
Cmc2i. The pattern weight factors proportional to the total number of counts 
per pattern were used. Quite similar results were found at varying the weight 
factor in the broad area of the 3D space of the weight factors. 

The use of multiple wavelengths allowed us to vary the contribution of cer- 
tain atoms to the diffraction patterns. Diffraction intensities around the angles 
of strongest peaks (29 ~ 30 — 40°) are modified by the anomalous scattering 
in such a way as if either Ho or Co were isomorphously replaced with La or 
Ti, respectively, in effect, as if the hypothetical compounds LaBaCo407 and 
HoBaTi407 were measured additionally to HoBaCo407. Fig. 5 shows the Ri- 
etveld plot of HoBaCo407 at the wavelength near Co K-edge. The inset shows 
the difference between two resonant patterns measured near the Co and Ho ab- 
sorption edges. These differential patterns consist of observed and theoretical 
profiles A/obs and A/theor at the profile conditions of the resonant pattern at 
Co K-edge. To plot the inset patterns the profiles Jobs and /thoor for A — 1.5421 
A were first transformed to the profile conditions of the profiles Jobs and /thcor 
for A = 1.6134 A and then subtracted from the latter. The difference between 
experimental and theoretical differential patterns is largest in a vicinity of the 
strongest peak given by the overlapped (023) and (113) refiections. In terms 
of both Pbn2i and Cmc2i structural models, the intensity of each of the (023) 
and (113) reflections should not vary with simultaneous changing the Co and Ho 
scattering power, as is expected from A/theor- However, it is observed in A/obs 
that the (023) reflection is enhanced at Co K edge, while the (113) reflection 
is suppressed. The whole difference profiles including the divergence between 
A/obs and A/thcor in the area of (023,113) reflections were observed to be nearly 
independent on the choice between the Pbn2i and Cmc2i structural models. 

Among the orthorhombic swedenborgites several compounds were refined 
previously with the S.G.'s Pbn2i and Cmc2i. The S.G. Cmc2i was found 
first to suit the structure of a Zn-based swedenborgitefisj . Later, the x-ray 
diffraction patterns of REBaCo407 for a number of RE were indexed with the 
orthorhombic cell of the Cmc2i symmetry [ij. We observed a few unindexed 
extra diffraction peaks, such as the reflection (105) in Fig. 2, which could be 
indexed within the S.G. Pbn2i suggested recently for YbBaCo407[3|. Since the 
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6Ba 
©Co 
•O 

Figure 6: Triple unit cell of HoBaCo407 refined with the symmetry group Cmc2i and viewed 
along [001]. The anisotropic thermal displacement factors refined for the atoms 031 and 033 
are represented by the surface of 50% probability. 



intensity of the extra reflections is extremely small, their occurrence has, in fact, 
a very small effect on the result of refinement. Therefore the R- factors are quite 
similar, and even better sometimes for Cmc2i than for Pbn2i. The results for 
both Pbn2i and Ctoc2i symmetries are shown in Table 1. Resulting structure 
in terms of Pbn2i is similar to that of YbBaCo4 06.95 at 150 kQ. The structure 
refined in terms of Cmc2i is shown in Fig. 6. 

In the higher-symmetry phase above Tg there are three nonequivalent oxygen 
sites. We call 01 the oxygens in the center of the star C04O13. Other oxygens 
at the bottom of K-layer in Fig. 1 are called 02 and those on top of K-layer are 
called 03. The atoms of cobalt and oxygen in the Tables 1, 2 and 3 are denoted, 
depending on the S.G., by single or double indices to simphfy the comparison 
between the low- and high-temperature phases. In the high-temperature phase, 
there occur just two inequivalent positions for Co labelled as Col, Co2. With 
lowering symmetry in the low-temperature phase, each of the Co2, 02 and 03 
splits into two (Ctoc2i) or three {Pbn2i) positions. Then the second index is 
added for these positions. 

Interatomic distances for the first coordination spheres in the low-temperature 
phase are shown in Table 2. Obtained for the S.G. Pbn2i{Cmc2i) average 
distances Col-0 and Co2-0 are 1.892(1.900) A and 1.918 (1.900) A, respec- 
tively. Average Co-0 distances for different Co sites in the kagome layer are 
(^)co2i-o = 1-921(1.958) A, (d)c„22-o = 1-907 A and (d)co23-o = 1-926(1.871) 
A. Some preference of the Col site for Co^+ ions can be suggested only for the 
S.G. Pbn2i, however, any charge stratification between K and T layers cannot 
induce orthorhombicity. Therefore the hypothesis of charge ordering is implau- 
sible within the Pbn2i model. In the Cmc2i model, the state with Co^+ placed 
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Figure 7: [001] projections of the crystal structure in high- (top panel) and low- (middle panel) 
phases in HoBaCo407 obtained in this work, and in YbBaCo407 (bottom panel, constructed 
for comparison according to the data of Ref . Q ) . 



into the smallest Co23 site is not a degenerated state because of double mul- 
tiplicity of this site compared to the multiplicities of the remaining Col and 
Co21 sites. Therefore, neither Pbn2i nor Cmc2i models fit the idea of charge 
ordering. 

We suggest another origin of the orthorhombicity, which is related to locking 
the rotational motion of tetrahedra in kagome layers. This motion is visualized 
with the use of anisotropic thermal factors in Figs. 6 and 7. Within the refine- 
ment of isotropic thermal factors the obtained i?iso were too unrealistic (> 4 A^) 
for several oxygen atoms in the structural models based on the S.G.'s Cmc2i 
(T < Tg) and PGsmc (T > Tg). Therefore, in these models, the anisotropic 
motion was allowed at fitting the thermal displacement factors of the 02 and 
03 atoms. This was not necessary for the Pbn2i model (T < Tg). One thus ob- 
serves at cooling through Tg locking the rotational motions in the 03 sites only 
for Cmc2i model, but in both 02 and 03 sites for Pbn2i model. In the high- 
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temperature phase PQ^mc, the correlated motion in the 02 and 03 sites can be 
thought as a rigid-body Ubration of the C0O4 tetrahedron. Such a rigid-body 
picture is not contemplable for the Cmc2i model, in which only the 031 and 
032 sites showed unusual thermal displacement factors. Related large deforma- 
tions of tetrahedra are unlikely. Therefore, through the rest of this article the 
Pbn2i model of the low-temperature phase will be considered in more detail. 



Table 2. Cation-oxygen bondlengths for all Ho, Ba and Co cations in 
HoBaCo407 at 300K. 



Bond 


Pbn2i 


Cmc2i 


Bond 


Pbn2i 


Cmc2i 


Bond 


Pbn2i 


Cmc2i 


H0-O2I 


2.291(13) 




Ba-023 


3.377(13) 


3.298(6) 


C02I-OI 


1.973(15) 


2.144(8) 


H0-O22 


2.152(12) 


2.190(8) 


Ba-031* 


2.819(10) 


3.162(5) 


C02I-O2I 


1.867(15) 




H0-O23* 


2.269(14) 


2.264(7) 


Ba-031 


3.494(10) 




C02I-O23' 


1.882(16) 


1.853(7) 


H0-O3I 


2.452(14) 


2.499(7) 


Ba-032 


2.940(14) 




Co21-O31 


1.963(17) 


1.98(1) 


H0-O32 


2.383(18) 




Ba-032 


3.387(14) 




C022-OI 


1.979(15) 




H0-O33* 


2.363(18) 


2.357(9) 


Ba-033* 


2.970(13) 


3.141(1) 


C022-O2I 


1.866(15) 










Ba-033'' 


3.353(13) 


3.19(1) 


C022-O22 


1.978(15) 




Ba-021 


2.742(11) 










Co22-032 


1.805(23) 




Ba-021 


2.534(12) 




Col-Ol 


1.897(8) 


1.910(9) 


C023-OI 


1.996(15) 


1.892(5) 


Ba-022 


2.969(12) 


2.785(8) 


C0I-O3I 


1.796(8) 


1.80(1) 


Co23- 022 


1.927(15) 


1.970(5) 


Ba-022 


3.313(13) 


3.498(8) 


C0I-O32 


2.011(17) 




Co23- 023 


1.977(15) 


1.883(8) 


Ba-023' 


2.900(13) 


2.979(6) 


C0I-O33* 


1.840(17) 


1.945(1) 


Co23- 033 


1.802(23) 


1.739(9) 



*The indicated bondlengths have double abundance in the S.G. C7Tlc2i due to the mul- 
tiplicity 8 for the positions of Co23, 023 and 033. 



The lattice parameters refined at T =380 K are a = 6.3046(2) A and 
c = 10.2604(4) A. The S.G.'s P31c and PQ^nic were examined and gave the 
quite similar parameters of fitting quality. In the Table 3, the atomic coordi- 
nates obtained within the S.G. PG^mc are shown. In a recent study of single 



crystals of HoBaCo407[lJ|, the arguments in favor of the latter S.G. were ob- 
tained from observing the equal intensities of symmetry-equivalent reflections. 
For example, the reflections (101) and (Oil) are symmetry-equivalent. Such 
distinguishable in single crystals reflections should have different intensities for 
P31c. In our Rietveld refinements with both S.G.'s P&^mc and P31c, we found 
while using i?iso, that the S.G. P31c permits to reach a better fit than the 
S.G. POsmc. However, with the S.G. PG^mc the equally good fitting qual- 
ity parameter (Psiagg = 5.65 ) was readily obtained through introducing the 
anisotropic factors of thermal displacements Py for 01 and 03 atoms. The 
refined anisotropic parameters for these atoms are shown as footnotes under 
Table 3. 

Table 3. Atomic coordinates for HoBaCo407 at 380 K. S.G. PQ^mc. 
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Atom 


X 


y 


z 


n. * 

-'--'ISO 


Ho 


2/3 


1/3 





0.4(3) 


Ba 


2/3 


1/3 


0.6258(7) 


1.5(3) 


Col 








0.569(3) 


0.4(2) 


Co2 


0.1705(7) 


-0.1705(7) 


0.805(2) 


0.4(2) 


01 








0.357(5) 


1.8(9) 


02 


0.498(3) 


-0.498(3) 


0.852(4) 




03 


0.837(3) 


-0.837(3) 


0.131(4) 





*Anisotropic thermal displacement factors were refined for Ol and 03. 
**{7il={/22=0.026{9), [733=0.07(1), C/i2=0.02(l), C/i3=-[/23=-0.005{5) A^. 
***[/ll = C/22=0.12{2), C/33=0.01(l), C/i2=0.11(2), J7i3=-{723=-0.002{4). 

Each of the Co sites in the kagome layer is larger than the Col site. 



In the high-temperature phase, we found the average distances Col-0 and 
Co2-0 of 1.96(2) A and 1.86(2) A, respectively. From low- to high-temperature 
phases, we observe, irrespective of the S.G. assumed for the low-temperature 
structure, that the average Co2-0 distances are shortened, while the average 
Col-0 distances are elongated. Especially, for the S.G. Phnlx each of the Co 
sites in the kagome layer is larger for T < Ts , but opposite is true for T > T^. 
The larger size for the Col site, {d)Q^i_Q > (1^)002-0 ^® ^^'^ argument against the 
larger valence of Col than that of Co2. The same trends were noticed by Huq et 
al 8], although these authors have refined the high-temperature structure with 
the S.G. P31c. We note, however, that the cobaltites REBaCo407, in which a 
larger valence was found for Col than for Co2 [§], did not show any structural 
phase transition. Judging from the relationship between ('i)cox-o ^^'^ ('^)co2-o 
such the transitionless cobaltites must be "in the low-temperature phase" with 
vanished by some reason structural distortion. The observed inversion of this re- 
lationship at heating through Tg has inspired our effort to search for the changes 
across Tg in the spectra of XAFS and Mossbauer. 

2.2. X-ray absorption spectra near Co K-edge 

The XAFS spectra of HoBaCo407 and Y'Qa.QoA-xFexO^ {x=Q and 0.1) are 
very similar (Fig. 8) The XANES begins with a pre-edge peak at E=7710 eV 
followed by the main edge and the " white line" maximum of absorption at 7730 
eV. These pre-edge peak and white line correspond to the electron transitions 
from Is orbital to the 3d and 4p bands respectively. Derivative spectroscopy 
using the second derivative of XAFS ensures the effective enhancement of reso- 
lution, which can be useful to separate two or more overlapping subbands. The 
pre-edge peak in tetrahedrally coordinated Co site is composed of the vacant 
orbitals situated in the e^ and t2g subbands. Even if they are unresolved in 
the pre-edge peak as wide as 5 eV, taking the second derivative of the XAFS 
spectra could be useful to resolve them as the separated minima of d'^fi/dE'^. 
In an octahedral Co site of Lai_a;Sra;Co03, for example, such a manipulation 
allows to resolve the e^ and t2g subbands separated by 2.5 eV 
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Figure 8: XANES spectra at the Co K-edge and the second XANES derivatives in the 
of the pre-edge peak for HoBaCo407 (1), YBaCo407 (2) and YBaC03.9Feo.1O7 (3). 
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Because our XAFS station was not equipped with a furnace, all measure- 
ments were done at 300 K. Instead of measuring the spectra of a REBaCo407 
across the phase transition, the spectra were measured for RE=Ho and Y with 
Ts=355 K and 280 K, respectively The only difference between the spectra 
for RE=Ho and Y is that the high-energy shoulder of d^^/dE^ near £'=7712 eV 
is better articulated in the low-temperature phase (RE=Ho). In a tetrahedral 
coordination, the Cg subband lies lower in energy than the t2g subband. There- 
fore, for the high-spin Co^+ {(f) the Cg orbitals are full and the t2g orbitals 
are half- filled. The predominant minimum of d^n/dE"^ should be ascribed to 
the t^g orbitals of the high-spin Co^+ associated with 75% of the total pre-edge 
peak area. Remaining 25% are due to Co^+. The K-edge Co'^"'" is expected to 
move to higher energy as a whole, including the pre-edge peak, by ca. 2 eV. 
In the high-spin (HS) state of the tetrahedrally coordinated Co'^"'", the splitting 
of the pre-edge peak is expected, since one half-filed orbital appears in the eg 
subzone, however, in the intermediate-spin (IS) state of such Co^"*", the pre-edge 
peak is unsplit again as it would come entirely from the transition Is — > 5d{t2g). 
Therefore, the feature near £=7712 eV must be better pronounced for the lower 
spin-state of Co^"*". Comparing Seebeck coefficient measured in YbBaCo407 
with theoretically predicted spin-dependent values Maignan et al.[^ found that 
both Co^+ and Co^+ are in HS states. Our spectra are consistent with the HS 
states of Co^+ and Co-^^, although a presence of some IS Co^+, especially in 
low-temperature phase, cannot be excluded. 

2.3. Mossbauer spectra 

Our Mossbauer spectra collected in the temperature range between ambient 
temperature and 133°C consist of two series. First, from one measurement to 
another, the temperature was ascending. In second series, the measurements 
were made for a sequence of temperatures at cooling the sample. All the ob- 
tained spectra are quite symmetric doublets (Fig. 9). Spectra treatments are in 
order, that would take into account the occurrence of four non-equivalent sites 
for the Fe dopants at T < Tg and two sites at T > Tg. However, up to date it 
was feasible only to fit each spectrum either with two-doublet or one-doublet 
(averaged) envelope. In the two-doublet spectra treatment, the following pa- 
rameters of the doublets were obtained at room temperature: S — 0.182 mm/s; 
£ = 0.389 mm/s and F = 0.342 mm/s for the majority (80%) subspectrum, and 
6 = 0.193 mm/s; e — 0.693 mm/s and F = 0.351 mm/s for the minority (20%) 
subspectrum. 

In the studied range of temperatures 297 K < T < 406 K, the average 
quadrupole splitting is varied by ~0.01 mm/s only. In such a narrow range, the 
determination of the behavior of the spectral parameters was feasible first of 
all due to high accuracy of our experimental setup. Second, a special analysis 
of the Mossbauer data treatments was performed. Namely, the accuracy of 
determination of spectral parameters in each spectrum was taken into account 
as described below. 

Parameters of the Mossbauer spectra were determined by a standard tech- 
nique via minimizing the value of ■ 
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Figure 9: A typical Mossbauer spectrum in HoBaC03.9Feo.1O7 (T= 300 K). 



N 2 
i=l 

i.e. finding the vector of parameters a of the theoretical spectrum Fth{i;a), 
such that = • In this work, for ah spectra we adopted the generahzed 
treatment, in which the parameter vector a was represented by four elements 
only: (1) average number of counts at infinite velocity, (2) chemical shift 5, (3) 
quadrupole splitting e; (4) linewidth of the lines of the doublet. Here i=l,...,N 
are the experimental points, Fex{i) is the number of counts in the i*'^ channel 
and (Ji is the evaluation of the dispersion in the experimental spectrum. 

During the spectra measurements in this work, the spectrometer operating 
time up to 1 month has been reached, however, the errors in determining the 
parameters from the minimization of in each spectrum (Eq.l) do not involve 
any drifts and/or uncontrolled instabilities related to a long-term experimental 
run. In order to estimate the relative contributions of the long-term sources 
of the errors the following procedure was adopted. First, we assumed that 
the chemical shift must depend on the temperature linearly as it does the so- 
called second-order-Dopplcr (SOD) shift, or "red shift". Therefore, we fitted 
27 obtained values of the chemical shift with the straight line. The errors of 
chemical shift obtained from non-linear regression of each of these 27 spectra 
(Eq.(l)) were used in this linear fit. As a result we obtained an estimation of 
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Figure 10; Temperature dependence of the chemical shift (5) and quadrupolar spHtting (e) 
in HoBaCog gFeo.oiOy. Open triangles shows the data obtained point by point sequentially 
at heating the sample; filled triangles shows the e-data obtained sequentially at cooling the 
sample. The straight line in upper panel is the linear fit to the experimental data: <5 = A + BT. 
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the squared deviation of the experimental chemical shift values (i.e., 5{n) as a 
component of the parameter vector a in Eq.(l)) from the linear temperature 
dependence hypothesis: 

X^({A^}) = E^^^^^^^4^^ (2) 

n=l " 

Here S=27, A=0.198 mm/s and B=-6.38-10-'* (mm/s) -K'^ are the parame- 
ters of linear temperature dependence of the chemical shift, obtained through 
minimizing the value of in the Eq.(2). The value of A refers to T=273 K. 

It was supposed then that the uncontrolled long-run instabilities of the spec- 
trometer can be taken into account through introducing a coefficient rj according 
to the following substitution A^^ — > v'^'^n- The value of rj can be estimated at 
the condition that the minimum of from Eq.(2) reaches its theoretical value 
when is replaced with ?y^A^: 

x'i{v'Al})^S-2 (3) 

Here S* — 2 is the number of degrees of freedom for the random variables 
in Eq.(2). The calculation performed according to this procedure has led to 
the result 77 = 1.6. Next, the errors in estimations of the chemical shifts and 
quadrupole splitting parameters were multiplied by the coefficient 77. The result 
of this operation is shown in Figs. (2) and (3). Finally, our analysis resulted in 
the average error of the quadrupole splitting not exceeding 0.0017 mm/s. The 
presented analysis is thus suggesting that the observed variation of e near Tg, 
although small, is the intrinsic and reproducible property of the material. 



3. Discussion 

The generalized procedure of analyzing the Mossbauer spectra allowed us to 
treat in the same way the spectra above and below Ts=355 K. More detailed 
treatments of our spectra, to be presented elsewhere, would take into account the 
occurrence of four non-equivalent sites for Co at T < Tg and two sites for Co at 
T > Tg. Recently, two subspectra were resolved at 4K in YBaCoQ2fiFeo.04O7.02 
owing to a difference of magnetic hyperfine fields AS^t — 3T[17|. The au- 
thors claim that the Mossbauer spectra cannot resolve the difference between 
the sites Co21, Co22 and Co23, however, it is not clear from their article 



17[, whether or not the orthorhombic distortion was present at all in their 



sample. No orthorhombicity down to 100 K was observed recently in single 



crystals of HoBaCo407[l^, however, the reason of such a suppression of the 
phase transition is yet not understood. When fitting our Mossbauer spectra of 
HoBaCo4Feo.i07 at 300 K (Fig. 9) with two doublets we observe the ratio of 
the areas of the doublets 4:1, with the 80% abundance of the narrower doublet. 
This is different from the ratio reported by Tsipis et al. (3:2) [17], although the 
radii of Y and Ho are similar, so that the distribution of Fe dopants over the 
Co sites is expected to coincide. However, the difference of distribution could 
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be related to the structural difference if the orthorhombicity in the samples of 
these authors[l7| is suppressed, similarly to the samples of Ref.|14|. Another 
reason for such a divergence may arise as a purely mathematical trick if the 
Mossbauer spectrometer instrumental linewidths and lineshapes are very dif- 
ferent. These arguments justify our choice in the present work of the method 
of spectra treatment, which simplify the comparison of obtained parameters 
between the samples with subtle structural differences. 

Both temperature regions above and below Ts=355 K in our HoBaCo407 
were characterized by very similar Mossbauer spectra. Also, the XAFS spectra 
of HoBaCo407 and YBaCo407 are very similar at 300 K despite Tg > 300 
K in HoBaCo407, but Ts< 300 K in YBaCo407[16(]. From these spectra we 
have shown that Co^"*" is in high-spin state. The same conclusion was drawn 
recently from the results of soft (Co ^2^3) x-ray absorption spectroscopy [Tsj. In 
addition, we have shown that the high-spin state is ubiquitous on both sides of 
Tg. Concerning the spin state of Co'^"'" the distortion of tetrahedral coordination 
would be associated to lifting the degeneracy either between double degenerated 
Cg states for S=2, or between triple-degenerated t2g states for S=l. For example, 
the S=l spin state correspond to a distortion which splits off a level for the hole 
location at the top of the band. Distortion of a tetrahedron via twisting a pair 
of vertices relative to other pair would produce such a splitting[2^ . In both low- 
and high-temperature phases, we cannot exclude the formation of intermediate- 
spin state for Co^+. Indeed, obtained distortions for the C0O4 tetrahedra are 



stronger than the distortions reported in Ref. [14| . 

Observed quadrupole splitting varies within 3% in the range ±AT/Ts==20%. 
The charge redistribution at Tg, seen not only via the change of the Co-0 
distances, but also via jumps of resistivity [3, looks after all very different 
than a charge ordering at a metal-insulator transition. Much larger changes of 
quadrupole splitting across Tg are usually expected when the charge ordering 
is driven by an orbital ordering. Huge jumps of e are typically associated with 
the orbital ordering. For example, in BiMno.95Feo.05O3, the jump-like ano maly 
As = 0.4 mm/s of quadrupolar splitting at Tqq — AlAK was observed[19j- 
In the manganites of the family REBaMn1.96Feo.04O6, fourfold increase of the 
quadrupole splitting was observed (at 300K) in the charge and orbitally or- 
dered compounds (RE=Y,Gd,Sm) compared to unordered ones (RE=La,Na,Pr) 
[iol . [ill , li^ . The orbital ordering coupled with charge ordering is therefore un- 
likely to drive the structural phase transition at Tg. 

Huq et al.[^ suggested the origin of the transition to be a "response to 
a markedly underbonded Ba^+ site in the high-temperature phase". The ion 
Ba^"*" occupies an anticubooctahedral site. In an ideal anticubooctahedral co- 
ordination, there are 12 equidistant ligands. The geometry of the anticubooc- 
tahedron close to ideal one was found in several studies on single-crystalline 
samples of REBaCo407 for RE=H oJl^ . LuQ, however, severely distorted in 
other structural studies on single-Qand poly-[l, [iJl crystalline samples. The 
Ba-0 distances of 3.14±0.01 A reported for nearly ideal anticubooctahedron by 
Juarez- Arellano et al. IJ] in HoBaCo407 and Kozeeva et al. 23[ in LuBaCo407 
seem to be too large for ordinary bondlengths of Ba^+ . No phase transition was 
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reported in such samples. In other studies, the distortion of the anticuboocta- 
hedron was reported to shorten the shortest Ba-0 bondlength down to 3.08 A 
in YBaCo407;Ti], 3.02 A in HoBaCo407 2i], 2.8 A in YbBaCo407[i]. The av- 
erage Ba-0 distances of 3.14±0.02 are nearly independent of RE, however, the 
distortions in much broader range were reported depending on the symmetry. 

In each of the structure models, PBamc, P31c, Cmc2i and Pbn2i, the distor- 
tion of the anticubooctahedron makes a number of Ba bondlengths increasing, 
and the same number of Ba bondlengths decreasing. The average Ba-0 distance 
remains unchanged. Namely, within the PQ^mc model, 6 oxygens of 12 remains 
at the ideal distances, 3 of them are allowed to move inwards and 3 other out- 
wards of the ideal anticubooctahedron. For the Cmc2i model, such a ratio 
changes from 3:6:3 to 5:2:5. The lower symmetries P31c and Pbn2i allow the 
distortion ratio 6:0:6, however, in YbBaCo407 Huq et alfs^ obtained within the 
P31c model the anticubooctahedron distortion scheme 3:6:3 similar to that of 
PG^mc model. This means the very small changes of bondlengths of 6 oxygens 
compared to the ideal anticubooctahedron. Therefore, the authors argued that 
the change of the distortion scheme of the anticubooctahedron from "3:6:3"-like 
above Tg to 6:0:6 below Tg is at the origin of the phase transition. This con- 
clusion is not contradictive to our viewpoint, the more especially as we fitted 
successively the high-temperature phase only with the PQ^mc model, which is 
intrinsically distorted according to the same 3:6:3 scheme. Although both the 
PG^mc and P31c models result in similar R-factors for the high-temperature 
phase, only the PG^mc model gave us plausible interatomic distances. Indeed, 
when trying to fit the high-temperature phase with the P31c model, we ob- 
tained the distortion of the anticubooctahedron much larger than that reported 
by Huq et al.Q. In addition, the unrealistic Co-0 distances of the order of 1.7 
A resulted, which obliged us to reject the P31c model. 

Thus the driving force of the transition to lower symmetry is likely to arise 
from too loose bonding of Ba^"*" in the symmetric anticubooctahedron. While 
the Ba^+ coordination polyhedron tends to collapse with lowering temperature, 
we do not observe that the coordination polyhedra of the Co ions change so 
dramatically. Displacive phase transitions in the tetrahedral networks are fre- 
quently driven by a soft optic phonon, which can propagate without distortions 
of tetrahedra. Such modes were called rigid- unit modes (RUMs) . 

The coordination polyhedron of Ba^+ is formed by 02 and 03 atoms. Each 
of these oxygens enters the first coordination sphere of two Co, two Ba^+ and 
one Ho^"*". On the other hand, the 01 atom enters to the first coordination 
sphere of four Co ions. Therefore, the deformation of the Ba^"*" coordination 
correspond to the motion of 02 and 03 atoms, which would be also represented 
in terms of kagome network as the motion of the of the 02 and 03 tetrahedral 
vertices around the 01 vertex. This is because each Co in the kagome layer 
(Co2 site) is coordinated by one 01, one 03 and two 02 atoms. 

Since the more covalent bonding is stronger within the framework of the 
C0O4 tetrahedra compared to the ionic bonding on the exterior of the frame- 
work, the distortion of the tetrahedra will cost more energy. The associated 
with such a distortion phonon modes would lie rather high on the frequency 
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Figure 11: The star-like fragment of the kagome layer showing the surrounding of the Ol 
" hinge" . 



scale. On the opposite, the RUMs that do not involve any significant stretching 
of the Co-0 bonds wiU have low frequencies. One of the low-energy RUMs ap- 
pears to be strongly temperature-dependent to become a soft mode driving the 
phase transition. This mode permits the intertetrahedral motions associated 
with bending the Co-O-Co links, where the central O atom can be viewed as a 
"hinge". The rotations of tetrahedra around the "hinges" are hindered either 
by the size of collapsed cavity around Ba^+ or by repulsion between oxygens 
getting too close to each other. Such hindered rotations named librations are 
commonly observed in tetrahedral networks [23|. In Fig.ll, the 0-" hinge" can 
be identified at the 01 atom (exactly at the center of the top panel in Fig. 7). 

From Fig. 7 it is evident that our low-temperature phase is quite similar to 
that reported for YbBaCo407[§|. It is also clear that the low-temperature phase 
can be obtained from the high-temperature one via quenching the positions of 
the 03 atoms at the endpoints of a counterclockwise rotation. Since the 03 site 
splits into three below Tg the rotation angles are different (in the range between 
5 and 10 deg.) for 031, 032 and 033. The rotation axis [001] is the same for 
031, 032 and 033, therefore, these oxygens keep their planar arrangement as 
in high-temperature phase (see the central plane at the bottom of T-layer in 
Fig.l). 

According to the RUMs concept, the axial rotation of the tetrahedral vertices 
located on top (in terms of Fig.l) of the kagome layers cannot occur without 
z-displacements of the 02 atoms at the bottom of these kagome layers. This 
is because only 25% of the oxygens at the bottom of kagome layers are located 
on the rotation axes. These 01 atoms constitute the rotation "hinges". Other 
75% oxygen sites in this plane are the 02 sites, which split into 021, 022 and 
023 in the low-temperature phase. Due to the rotation of the tetrahedral "rigid 
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units" the 021 and 022 atoms get out of the plane, in which they were in the 
high-temperature phase. 

The shift of the 021 site along z-axis is as large as 0.52 A. This is consistent 
with the largest rotation angle (w 10°) for the 031 atom, which belongs to the 
tetrahedra of Co21 (in Fig. 7 two of these tetrahedra form a vertical "sandglass" 
figure extended along 6-axis around the central 01 atom). Smaller rotation 
angles are shown in Fig. 7 for the 032 and 033 atoms, which belong to the 
tetrahedra of Co22 and Co23, respectively. This is in agreement with smaller 
out-of-plane 2-shift of the 022 atoms (0.14 A), which belong to both Co22 and 
Co23 tetrahedra. The third site 023 remains at the same z-level as in high- 
temperature phase. This site belongs to Co21 and Co23 tetrahedra. 

Thus, the tetrahedra faces, which were parallel in high-temperature phase, 
become slightly tilted relative to horizontal plane. The measure of tilts for the 
bottom faces of Co21 and Co23 tetrahedra are the z-shifts of corresponding 
vertices, i.e., 0.52 A and 0.14 A, for Co21 and Co23 tetrahedra, respectively. 
The degree of tilt for the Co22 tetrahedra, which is formed by the 021 and 
022 vertices is measured by the difference of their z-shifts, is 0.52 — 0.14 = 
0.38 A. This relationship directly follows from the symmetry Pbn2i. Because 
of such a relationship between the values of z-shifts for 021, 022 and 023 
atoms the spread between the tilts of three tetrahedra surrounding the 01- 
hinge becomes unavoidable. Quite similar spread was obtained by Huq et al.Q. 
If the tetrahedra are indeed behaving as rigid-units, the related spread would in 
turn appear in the angles of rotation of the three 03 atoms around the 01 — z 
axis. Comparing the results of our refinement with the data of Ref. 0], we show 
in Fig. 10 that the spread is quite similar for HoBaCo407 and YbBaCo407. 

In high-temperature phase, where all three 03 atoms around the 01 — z axis 
are equivalent, they librate with equally large amplitudes as shown at the upper 
panel of Fig. 7. The thermal displacements of 03 atoms have no component 
along z, however, the z— component is the main component for 02 atoms (see 
the footnote of Table 3). However, because the rigidity of tetrahedra put the 
constraint on the z— shifts of three 02 atoms, the rotations cannot be condensed 
below Ts within the high-symmetry model PG^nic, in which all of the 02 sites 
are equivalent. Therefore, the rigid-body libration of tetrahedra lowers the 
symmetry, so that the orthorhombicity directly follows from combining the RUM 
condensation with translation symmetry. 

In tetrahedral networks, various constraints on possible models of deforma- 
tion of these networks were described previously by Dove et al[2l|. In the 
anion-close-packed systems, the structure becomes highly susceptible to the 
libration-like motions when the proportion of " bridgelike" two-coordinated an- 
ions is large. For example, in the red mercury Hgl2, where all the anions are 
bridgelike, the coupled librations of IIgl4 tetrahedra are unconstrainted and 
the thermal displacements of anions are strongly anisotropic [2^. This is not 
the case for the built of tetrahedra and octahedra structure of spinel M3O4, in 
which all the oxygens are four-coordinated. A detailed comparizon of the struc- 
ture of REBaCo407 with spinels, hexaferrites, perovskites and with YBaFe407 
ferrite was recently made by Caignaert et al[30|. In these structures, the fraction 
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of bridgelike oxygen increases with increasing the fraction of large-size cations. 
Like in the red mercury, in the perovskites, all the oxygens are bridgelike and 
the perovskites are reach by their famous soft modes. 

Phase transitions known for the tetrahedral and octahedral networks are 
classified into displacive transitions and the transitions of the order-disorder 
type. In a simple model of double- well potential[l^, the type of transition 
is determined by the relationship between energy of harmonic forces between 
neighboring atoms A:bTs) and the height of the barrier Vq between two min- 
ima of the potential well. The degree of displaciveness for a transition is given 
by the parameter = fce^s/Vb- We are in the displacive hmit when :s> 1 
and in the order-disorder limit when s 3> 1. Interestingly, an intermediate type 
of transition (s ~ 1) was attributed to freezing the librations of large amplitude, 
such as librations of CO3 groups to angle ~ 30° in calcite[lo|. 

In the displacive limit, there exist a well established relationship between 
the transition temperature and the angle of rotation of tetrahedra at K: 

ksTs = Kifl (4) 



Markina et al.Q and Juarez- Arellano et al. 14 1 have reported nearly linear rela- 
tionship between the value of Tg in REBaCo407 and the size of ionic radii of RE. 
The temperature Tg increases twice when the size of RE increases by 4%. The 
radii of Yb and Ho differ by 3%. However, in Fig. 7 we observe that the values of 
(fio are very similar for RE=Yb and Ho. The force constant K in the Eq.(4) is 
the stiffness of the rigid-body unit. Therefore, if we are in the displacive limit, 
where the Eq. (1) is valid, the linear relationship between Tg and K may signify 
the increase in the strength of the Co-0 bonds with increasing the ionic radius 
of RE. This phenomenon can be regarded as the manifestation of the inductive 
efFect[3l| in the kagome network: Co-0 bonding is reinforced at the expense of 
RE-0 bonding. Indeed, all the 02 and 03 sites (85.7% of all oxygen sites) are 
coordinated by two Co, two Ba, and one RE. The inductive effect[31] related 
to decreased covalence of RE-0 bonding with increasing the size of RE would 
concern these 02 and 03 sites. 

Only remaining 14.3% oxygen sites (01-sites) are not bound to RE, as they 
are coordinated by four Co. With respect to the covalent Co-0 bonding, a 
bridge-like (02 and 03) and tetrahedral (01) oxygen can be distinguished. The 
latter shows more ionic bonding and longer Co-0 bonds. 

All Co ions are, in their turn, containing in their first coordination sphere 
three bridgelike oxygens and one tetrahedral oxygen. Two types of ligands for 
each Co presents the plausible argument to explain rather narrow lines of the 
symmetric doublet observed in the Mossbauer spectra. In spite of the occurrence 
of four different sites for Co in the low-temperature phase, the spectrum shown 
in Fig. 9 evidences the narrow distribution of EEC over these sites. If we 
attribute the origin of EEC to strong difference between the bridgelike and 
tetrahedral oxygen ligation, then the EEC can be explained in terms of point 
charge model for the first coordination sphere. Then due to the symmetry in 
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high-temperature phase the EFG on each of the Co sites should be oriented 
along the longest Co-01 bonds for both Col and Co2 sites. 

The tetrahedra are non-ideal rigid bodies, therefore, in low-temperature 
phase each of three tetrahedra in kagome layer is distorted in its own manner. 
However, we observe only the change of e within 3%. This is an important ob- 
servation, which would let us to estimate whether or not the partial contribution 
of each ligand to EFG depend on the bond angles. Generally, upon variation 
of bond angles the partial contribution to EFG may change because fractional 
p— character of the bonds changes [s^ . However, the families of compounds ex- 
ist, in which the assumption of angle- independent partial contributions is a good 
approximation. One of such families includes, for example, brownmillerites and 
high-Tc superconductors with "1212" structure [11]. The high-Tc cuprates, so 
different at first sight from brownmillerites, show very similar ligation for iron 
located at the midway between the cuprate planes 33, H, 3^. 



The quadrupole splitting is related to the EFG tensor main component 



and EFG asymmetry rj = {Vxx — Vyy)/Vzz viaj35l|: 



2E^ 

Here = 14.4125 keV is the transition energy, and Q is the quadrupole mo- 
ment of the excited state of the ^^Fe nuclei. From the electric field gradient 
Yzz in units of V/m^ the quadrupole splitting in mm/s can be obtained using 
the factor i • eQcE''^ = 0.1664 • lO-^imm-s-^V-im^ for the Q value of 0.16 
barn. The Eq.(5) shows that the changes in e signifies either variation of Vzz or 
variation of rj or both. 

The behavior of r] is well predictable in a vicinity of a second-order phase 
transition. Displacive phase transitions are usually close to the continuous phase 
transitions, which are described though rigorously with renormalization group 
theory as first-order-like, in practice exhibit the size of first-order discontinu- 
ities so small as to be virtually unobservablejlOj . Through the continuous phase 
transitions an empiric rule was established[36|, which states that all three com- 
ponents of EFG tensor depends linearly on a single control parameter. In other 
words, the changes of principal component and asymmetry of EFG tensor are 
correlated and both Vzz and 77 vary continuously through such a transition. 
Symmetry breaking obviously underlies this rule in a vicinity of axiality (77 = 0) 
or antiaxiality {jj = 1). 

The non-trivial correlation between Vzz and r/ is illustrated by the linear 
trajectories of the system in the so-called Czjzek-plot, in which — 2V^x is plot- 
ted versus | Kz | (/3 -I- r7//3)[3l, Is^]. These axes are chosen to make use the 
intrinsic properties of the EFG tensor, which is symmetric and traceless. There- 
fore at the condition of antiaxiality (at horizontal axis in Czjzek plot, Fig. 12a) , 
where Vzz changes its sign by definition, the system pass the line rj ~ 1 without 
refraction. On the other hand, at 77 = the system reflects from the borders 
of Czjzek plot with billiard kinematics. These properties are satisfied always 
when each component of the EFG tensor depends linearly on a single control 



24 



E 
> 

CM 

o 



> 

CM 




1 2 3 4 5 6 

2/V3 I 2V +V I (10'' V/m') 



1 
1 







>' 
II 
^ 


-0 



> 
> 



1 

OH 

9 
8 

7-\ 
6 

5H 

4- 

3A 

2 

1 

OH 
1 



- 771/12 



Un/24 




5ji/8 to 
-7Tt/24 



-5;t/12 



20 40 60 80 100 120 140 160 



T(°C) 



Figure 12: (a) Czjzek-plot of Y = Vzz{l - v) = -SVxx versus X =| | (</3 + r;/^3) 
= 2/v^3 I 2K!z + Vlra; I with the EFG tensor components ordered according to | Vxx |<| 
Vyy \<\ Vzz I- Lines originating from the origin are rj = const, lines. The angles of inclination 
of the straight-line trajectories with respect to abscissa axis are indicated. Arrows show the 
direction of increasing temperature, (b) Corresponding dependences of the EFG asymmetry 
TjiT). The frames mark the 77(T)— curves for two trajectories whose EFG components are 
shown in detail in Fig. 13. 
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Figure 13: Variation of the EFG tensor components for two straight-line billiard trajectories 
in Czjzek-plot with the direction of increasing temperature inclined at the angle to abscissa 
axis of 57r/8 and reflected to —7it/2A (a) and 137r/24 (b). Corresponding dependences r]{T) 
are marked by frames in Fig. 12 (b). 
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parameter p[3a]- This parameter is in turn a function of temperature p = F{T). 

In Fig. 12(b), we plotted r] versus T assuming several linear paths of the 
system in Chjzek-plot shown in Fig. 12(a). First of all, the path at the angle 
of 7r/2 to horizontal axis crosses the line 77 = 1 at the same temperature where 
e(T) shows the minimum in Fig. 10. Therefore, rj{T) curve culminates at 
T = 70°C in Fig. 12(b). Also, the trajectory, which starts at the angle of 27r/3 
near 77 = 0.4, then goes to the upper border of the Chjzek-plot, reflects from 
the 77 = line backwards to the angle of -7r/3 and returns back to 77 = 0.4, 
produces the curve 77(T) with the minimum at T = 70°C. The path tilted to the 
midway angle with the directions towards (-57r/12) and (77r/12) produces two 
77(T) curves with no extremum, but with sigmoidal shape having the inflexion 
again at T = 70° C. These two curves correspond to two opposite directions of 
increasing temperature along the same line indicated by arrows on Czjzek-plot 
(Fig.l2,a). 

Next, we investigate the trajectory tilted to the angles intermediate between 
77r/12 and 27r/3. As in the case of 27r/3, such a trajectory is reflected from the 
border of Czjzek-plot, but now not at the right angle. The path starting at the 
angle of 177r/24 reflects to the angle of — 37r/8. Moving in the opposite direction 
along the same trajectory we start from short segment at the angle of Stt/S and 
reflect to -77r/24. In both these cases, the reflection point does not coincide with 
the minimum of the £(T) curve in Fig. 10. Because our experimental data for 
e(r) are more or less symmetric relative to the minimum at T = 70°C, two 
possible trajectories for these angles generate two curves for r}{T). One of them 
shows minimum at T = 82 °C and another shows minimum at T = 60 °C. 

In the same way, two trajectories in the Czjzek-plot and two 77(r) curves can 
be generated for the angles of -ll7r/24 and 137r/24 along the line intermediate 
between 77r/12 and 7r/2. These two paths cross the horizontal axis (r/ = 1 line) 
at different temperatures. Again the generated rj{T) curves culminate at T = 82 
°C and at T = 60 °C. 

Since we know now both values of 77(T) and | V^ziT) \ along each path in 
Czjzek plot, we can derive each component of the EFG tensor. In the displacive 
limit, the behavior of the 77(T) in a vicinity of Tg is expected to follow closely the 
behavior of order parameter around the point where symmetry breaks 36| . Two 
of the obtained 77(T) curves culminate indeed at Tg = 82 °C. Therefore, if we 
are indeed not far from the displacive limit, the EFG components T4z, Vyy and 
Vxx should behave as shown in Fig. 13. In the first case (57r/8 to -77r/24), the 
EFG is nearly axial above Tg — 82 °C, but quickly approaches to antiaxilality 
as the hexagonal symmetry breaks down. The non-principal EFG components 
exchange their direction at Tg. The second case (137r/24) is opposite. In the 
second case, the principal EFG component switch its direction at Tg. The plot in 
Fig. 13(a) corresponds to the appearance of the additional source of EFG below 
Tg, which induces second large EFG component. On the opposite, the plot in 
Fig. 13(b) correspond to the case when two large components of EFG exist above 
Tg and one of them is extinguished at the expense of the displacements showing 
up with the breakdown of symmetry. 

Only the first case fits the high-temperature axiality which should correspond 
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to OUT classification of the oxygen ligands as bridgelikc and tetrahedral ones. In 
this case the z-axis of the EFG in its internal coordinate system must be directed 
from each Co ion towards ionic tetrahedral oxygen 01 (in the radial direction 
of Fig. 11). In this picture, the small asymmetry of EFG above Tg arises from 
the structural difference between bridgelike oxygens 02 and 03. Because of 
non-ideal rigidity of tetrahedra below Tg the displacements of the 03 atoms 
arc larger than the displacements of the 02 atoms. The change in the bond 
angles and the associated change in the fractional p— character of these oxygen 
bonds result in the additional component of EFG, which becomes comparable 
in magnitude with the main EFG component. 

4. Concluding Remarks 

We have studied the temperature dependence of structure parameters and 
measured the hyperfine splitting e(T) of ^''Fe in HoBaC03.9Feo.1O7 in the range 
±AT/Ts=:20% around the structural phase transition. Although the changes of 
£(Ts± AT) within only 3% were observed, we showed that the non-principal EFG 
components {Vxx and Vyy) vary much more dramatically in a vicinity of Tg = 355 
K in agreement with the local structure changes at the symmetry-breaking tran- 
sition. Although the variations of each of three individual EFG components Vxx, 
Vyy and Vzz are smooth similarly to the temperature dependences of the lattice 
parameters a and 6/^/3, we have demonstrated the nonmonotonic behavior of 
the combination r]{T) = {Vxx — Vyy)/Vzz, which experiences a dip-like (or cusp- 
like) anomaly at Tg. In some vicinity of Tg the temperature dependence of the 
EFG asymmetry ??(T) is that of the order parameter as expected for a displacive 
phase transition. 

Rietveld refinements of the structure on both sides of Tg allowed identifying 
the main component of thermal motion which freezes as the hexagonal symmetry 
breaks down. The origin of the transition is ascribed to the condensation of the 
libration phonon mode associated with the rigid-body rotational movements of 
starlike tetrahedral units, the building blocks of kagome network. 
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7. Figure Captions 

Fig. 1. Perspective view of the quadruple unit cell in the low-temperature 
phase of ReBaCo407. Kagome and triangular layers of C0O4 tetrahedra are 
marked by "K" and "T" . The quadruple unit cell of HoBaCo407 is shown with 
atomic coordinates refined in S.G. Pbn2i. 

Fig. 2. Three wavelengths used in the diffraction experiments shown on the 
plot of theoretical atomic scattering factor corrections // and /" for Co and 
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Ho (a); characteristic area of the high-resolution synchrotron x-ray diffraction 
patterns (Rietveld plots) at A = 1.5421 A in the high-temperature (b) and 
low-temperature (c) phases of HoBaCo407. Structure of the low-temperature 
phase was refined at T = 300K using either S.G. Pbn2x or Cmc2i, one of 
which {Pbn2x) is shown in (c). Theoretical positions of the permitted reflections 
are shown in (c) for both Phn2x and Ctoc2i groups. Indicated by the arrow 
reflection (105) is among very weak peaks extinguished for S.G. Ctoc2i, but 
allowed for S.G. P6n2i. 

Fig. 3. The curves of differential thermal analysis in HoBaCo407 measured 
upon heating and cooling at a rate of 3 deg./min. 

Fig. 4. Evolution of the lattice cell parameters a, h/^'i and c with temper- 
ature near the structural phase transition. 

Fig. 5. The conventional (main panel) and differential (inset) Rietveld plots 
for the synchrotron x-ray diffraction data collected at T=300 K (below the phase 
transition) in HoBaCo407. Main panel: observed, calculated and difference 
intensities near Co K-edge (A — 1.6134 A). Inset: the differences /obs(CoK) — 
/obs(HoLiii) and /thcor(CoK)- /thoor(HoLiii), where the profile /thcor(HoLiii) 
was transformed to the /thcor(CoK) profile conditions. 

Fig. 6. Triple unit cell of HoBaCo407 refined with the symmetry group 
Ctoc2i and viewed along [001]. The anisotropic thermal displacement factors 
refined for the atoms 031 and 033 are represented by the surface of 50% prob- 
ability. 

Fig. 7. [001] projections of the crystal structure in high- (top panel) 
and low- (middle panel) phases in HoBaCo407 obtained in this work, and in 
YbBaCo407 (bottom panel, constructed for comparison according to the data 
ofRef. B). 

Fig. 8. XANES spectra at the Co K-edge and the second XANES derivatives 
in the region of the pre-edge peak for IIoBaCo407 (1), YBaCo407 (2) and 
YBaC03.9Feo.1O7 (3). 

Fig. 9. A typical Mossbauer spectrum in IIoBaC03.gFeo.1O7 (r= 300 K). 

Fig. 10. Temperature dependence of the chemical shift [5) and quadrupolar 
splitting (e) in IIoBaCo3 gFeo.oi07. Open triangles shows the data obtained 
point by point sequentially at heating the sample; filled triangles shows the 
data obtained sequentially at cooling the sample. The straight line in upper 
panel is the linear fit to the experimental data 5 ^ A-{- BT. 

Fig. 11. The star-like fragment of the kagome layer showing the surrounding 
of the 01 "hinge". 

Fig. 12. (a) Czjzek-plot of F = 14^(1 - r/) = -2t4x versus X =| t4z | 
(/3 + r]/'/3) = 2//3 I 2Vzz + Vxx \ with the EFG tensor components ordered 
according to | V^^ [<[ Vyy [<[ t4z [• Lines originating from the origin are 
77 = const, lines. The angles of inclination of the straight-line trajectories with 
respect to abscissa axis are indicated. Arrows show the direction of increasing 
temperature, (b) Corresponding dependences of the EFG asymmetry 77 (T). The 
frames mark the 77(T)— curves for two trajectories whose EFG components are 
shown in detail in Fig. 13. 

Fig. 13. Variation of the EFG tensor components for two straight-line 
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billiard trajectories in Czjzek-plot with the direction of increasing temperature 
inclined at the angle to abscissa axis of 57r/8 and reflected to — 77r/24 (a) and 
137r/24 (b). Corresponding dependences ri{T) are marked by frames in Fig. 12 
(b). 
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